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^ ■ Abstract 

QhI We discuss to what extent future precision measurements of neutrino mixing 

Q^' observables will influence the information we can draw from a measurement of (or an 

1^ . improved limit on) neutrinoless double beta decay. Whereas the Am^ corresponding 

to solar and atmospheric neutrino oscillations are expected to be known with good 
precision, the parameter 9i2 will govern large part of the uncertainty. We focus in 
^ ■ particular on the possibility of distinguishing the neutrino mass hierarchies and on 

setting a limit on the neutrino mass. We give the largest allowed values of the neutrino 
masses which allow to distinguish the normal from the inverted hierarchy. All aspects 
are discussed as a function of the uncertainty stemming from the involved nuclear 
matrix elements. The implications of a vanishing, or extremely small, effective mass 
are also investigated. By giving a large list of possible neutrino mass matrices and 
their predictions for the observables, we finally explore how a measurement of (or an 
improved limit on) neutrinoless double beta decay can help to identify the neutrino 
mass matrix if more precise values of the relevant parameters are known. 
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1 Introduction 



The last few years have seen tremendous progress in our understanding of the neutrino 
sector. The solar neutrino deficit is now known to be certainly due to neutrino flavor 
oscillations [UlSElllllSlEl with the best-fit oscillation parameters given as Am^ = Aml-^ = 
8.0 X 10~^ eV^ and sin^ 9q = sin^ = 0.31 The atmospheric neutrino deficit can also 
be ascribed to flavor oscillations with a good degree of confidence [3 |S] with the best-fit 
oscillation parameters as Am^tm = |Am|^| = 2.1 x 10~^ eV^ and sin^ 2^atm = sin^26'23 = 1 
[Zj. The upper limit on the third mixing angle ^13 is mainly determined by the reactor 
neutrino data [H], which when combined with the information obtained from the solar and 
atmospheric neutrino experiments gives a bound of sin^ ^13 < 0.044 at 3a piU lllj. 

Despite all these impressive achievements, there still remains a lot to be learned. The 
arguably most fundamental question, namely whether neutrinos are Dirac or Majorana 
particles, remains still to be answered. If neutrinos are Majorana particles, we have nine 
physical parameters describing the 3x3 light neutrino mass matrix and determining them 
is one of the ultimate goals of neutrino physics. With the two Am^ and two mixing angles 
- 612 and ^23 — known, there still remain five parameters to be tracked down. To be more 
precise, we still lack knowledge about the following points: 

• Is there CP violation in the lepton sector in analogy with that in the quark sector? 

• Is the atmospheric neutrino mixing angle ^23 exactly maximal? If not, does it lie 
above or below it/ 47 

• Is the third mixing angle ^13 exactly zero? 

• How small is the absolute neutrino mass scale? 

• What is the ordering of the neutrino masses, i.e., what is the neutrino mass hierarchy? 

Answers to all or some of these question will certainly lead to better understanding of the 
underlying theory which gives rise to neutrino masses^. 

Neutrinoless double beta decay (Oz//3/3) is expected to be of crucial importance in an- 
swering some of the questions raised above. At the same time, answers to some the ques- 
tions obtained from elsewhere will help in interpreting a positive or negative signal in the 
next generation Oz//3/3 experiments. Goal of these experiments is the observation of the 
process 

{A,Z) {A,Z ~2) + 2e- . 

rarely mentioned point is that for Majorana neutrinos one can, at least in principle, fully determine 
the complete mass matrix because it is symmetric. If the neutrinos are Dirac particles like quarks and 
charged leptons, their mass matrix is in general not symmetric and therefore parametrized by fifteen 
free parameters. However, since we only see left-handed weak currents, the "right-handed" parts of the 
neutrino parameters are out of reach and hence not observable. 
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The effective mass which will be extracted or bounded in a Oi//3/3 experiment is given by 
the following coherent sum: 



{m)= Y^^^Ul , (1) 

i 

where rrii is the mass of the i*'' neutrino mass state, the sum is over all the light neutrino 
mass states and Uei are the matrix elements of the Pontecorvo-Maki-Nakagawa-Sakata 
(PMNS) neutrino mixing matrix jT^. We parametrize it as 



U 



( C12C13 S12C13 Si3 

— ■512C23 — Ci2S23Si3e*'^ C12C23 — Si2S23Si3e*'^ ■523Ci3e"^ 
V ■S12S23 — Ci2C23Si3e"^ — C12S23 — Si2C23Si3e*'^ £230136*'' / 



diag(l,e^^e^^), (2) 



where we have used the usual notations Cij = cos 6ij, Sij = sin 6ij, 6 is the Dirac CP- 



violation phase, a and /? are the two Majorana CP-violation phases Thus we can 

define (m) in terms of the oscillation parameters, the Majorana phases and the neutrino 
mass scale. This means that (m) depends on 7 out of 9 parameters contained in the 
neutrino mass matrix. In particular, the effective mass is a function of all unknowns of 
neutrino physics except for the Dirac phase and 6'23. This means that a measurement of - 
or a stringent limit on - OujSjS will give us some information on the unknown parameters. 
Combined with complementary information from other independent types of experiments, 
this would reinforce our understanding of the neutrino sector. 

There is a large list of available analyzes focusing on the connection of the known and 
unknown neutrino parameters with the effective mass fSl ^1 ^1 QHl 120], for a recent 
review of the theoretical situation see [21]. In this work we want to focus on the mass 
ordering of the neutrinos, on the neutrino mass scale and on the implications of a very 
small or vanishing effective mass. The effective mass to be extracted from neutrinoless 
double beta decay depends crucially on the neutrino mass spectrum. Of special interest 
are the following three extreme cases: 



normal hierarchy (NH): |m3| ~ \/ /S.m\ ^ |m2| — A/AmQ ^ |mi| , (3) 



inverted hierarchy (IH): |m2| — \mi\ ~ y^AmJ ^ , (4) 



quasi-degeneracy (QD): mo = \'m^\ — \m2\ — \mi\ ^ ^J Am\ . (5) 



The order of magnitude of the effective mass in those spectra is ^y ArriQ, \/ Am^ and mo, 
respectively. 

The best current limit on the effective mass is given by measurements of ^^Ge established 
by the Heidelberg-Moscow collaboration [22] 

(m)<0.35^eV, (6) 



3 



where z = 0{1) indicates that there is an uncertainty stemming from the nuclear physics 
involved in calculating the decay width of Ouf^p. Similar results were obtained by the 
IGEX collaboration |23I. Several new experiments are currently running, under construc- 
tion or in the planing phase. The NEMOS ^24j and CUORICINO [22] experiments are 
already taking data and reach sensitivities near the current limits. The next generation of 
experiments, with projects such as CUORE [21, MAJORANA |27], GERDA [28j, EXO 
[2ni, MOON jSO], COBRA jSH, XMASS, DCBA 32 , CANDLES CAMEO P], (for 
a review see jHH]) will probe values of (m) one order of magnitude below the limit from 
Eq. (jni). Thus we expect (m) to be probed down to ~ ■>/ Am\ ~ 0.04 eV and it would 
be pertinent to ask if such a measurement could help us learn about the neutrino hierarchy^. 

Maximal mixing in the solar neutrino sector is now disfavored at almost 6a and this 
non-maximality of 612 allows for the possibility to distinguish the NH from the IH case. 
This however depends on the value of the smallest neutrino mass state. We investigate this 
issue and give the value of the smallest neutrino mass state which allows to distinguish NH 
from IH. Another aspect, concerning the QD spectrum, is the size of the absolute neutrino 
mass scale. We give and analyze a simple formula for a limit on mo in terms of 612, 613 and 
the limit (or value) of (m). All these issues depend crucially on the uncertainties involved. 
In addition to the experimental errors, we would have to contend with the theoretical un- 
certainties coming from unknown nuclear matrix elements and uncertainties in the values 
of the oscillation parameters. As far the oscillation parameters are concerned, the uncer- 
tainties are expected to be reduced sharply by the next generation experiments. This will 
greatly reduce the expected range of (m) for a given mass spectrum. We give a detailed 
review of the improvements expected in each of the oscillation parameters and identify the 
parameters which would still lead to maximum uncertainty in the predicted value of (m) 
for a given mass spectrum. We take the uncertainties coming from our lack of knowledge 
of the nuclear matrix elements and explore the chances of determining the neutrino mass 
hierarchy from a positive Oi^Pf] signal in the future. Furthermore, we analyze the conse- 
quences of a very small or even vanishing (m) if neutrinos are indeed Majorana particles. 
Such a small (m) will influence the allowed range of values of the mixing parameters 612, 
6'i3, the Majorana phases and the smallest neutrino mass. We compare these implied range 
of parameter values with their current limits. Finally, we perform a thorough scan of the 
phenomenologically viable neutrino mass models and give a list of possible neutrino mass 
matrices out of which the true one might be identified when future precision measurements 
of neutrino oscillation observables are combined with a measurement of or limit on (m). 

We begin by giving in Section |21 a detailed report of the existing range of values for the 
oscillation parameters. We list down the expected improvements in each of the parameters. 
In Section El we probe the chances of distinguishing the IH from the NH from a future posi- 
tive signal for neutrinoless double beta decay. We take into account the current and future 

^Of course those experiments aim also to put the controversial evidence of part of the Heidelberg- 
Moscow collaboration to the test. 
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expected uncertainties on the oscillation parameters and include the role of the uncertainty 
in the nuclear matrix elements. In Section |3] we consider QD as the true mass spectrum of 
the neutrinos and study what limits one could set on the absolute/common neutrino mass 
scale mo from a measured value of (m) and the current and future uncertainties on the 
oscillation parameters. Sectional probes the situation one would have if we do not see any 
signal in the next set of Oz//5/5 experiments and yet believe that neutrinos are Majorana 
particles. In Section IHl we then list viable neutrino mass matrices and discuss how the 
number can drastically be reduced when measurements of (m) and the other oscillation 
parameters have been performed. We finally conclude in Section 

2 Past, Present and Future data 

In this Section we shall present the status of the current global neutrino data and the 
prospects for its future improvement. First we give the current situation. 

At the 3a level, our current knowledge of the solar parameters within 3a is limited to 

lanni 

7.0 X 10-5 eV^ < Am| < 9.3 x lO'^ eV^ , (7) 
0.24 < sin^ 012 < 0.41 . (8) 

The atmospheric mass squared difference and mixing angle at 3a are known within (71 

1.3 X 10-3 eV^ < Ami < 4-2 x lO'^ eV^ , (9) 
0.33 < sin^ 023 < 0.66 , (10) 

while the mixing angle 6*13 at 3a is restricted to lie below the value UHl HH ISI] 

sin^ < 0.044 . (11) 

A very useful parameter related to the mass hierarchy of the neutrinos is the ratio of 
the solar and atmospheric mass squared differences, 



R = — ^ , with 0.017 <R< 0.072 at 3a . (12) 



Regarding the absolute value of the neutrino masses, several approaches exist. The most 
model-independent one is certainly the direct search for kinematical effects in the energy 
spectra of beta-decays. The Mainz [3H] and Troitsk [32] experiments gave upper limits 
on the electron neutrino mass of 2.3 eV at 95 % CL. Cosmological observations imply 
typically more stringent limits, which however considerably depend on the data set and 
the priors used in the analysis. Combining the cosmic microwave radiation measurements 
by the WMAP satellite [101 with data on the large scale structure of the Universe and 
other data sets, gives limits on the sum of neutrino masses between 0.4 and 2 eV, see ^2] 
for a recent update of the situation. 
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Finally, there is currently no information on any of the three possible CP phases, which 
take on values between zero and 2tt. 

Let us now discuss the current uncertainties in the values of the parameters and their 
expected improvement. We can see from Eqs. ((Zj) and © that we still have 53% (14%) 
uncertainty^ on the value of Amg^^ (Am2i) at 3a. These ranges of allowed values are 
expected to reduce remarkably in the future: the uncertainty in Amgi is expected to reduce 
to 7% at the 3a level from future measurement of reactor antineutrinos in KamLAND 
inZlE] and the value of Am|^ is expected to be known with ~ 5% at the 3a level from the 
next generation accelerator experiments involving superbeams in the next ten years |42j . 

The uncertainty in the knowledge of sin^ 623 is at present 33% at 3a. This is expected to 
reduce to ~ 20% in the next ten years after the results from the next generation superbeam 
experiments become available It is clear that the mixing angle 623 plays no role for 
the study of neutrinoless double beta decay. Nevertheless, we shall see later on in Section 
|H]that in order to identify the neutrino mass matrix it can be rather important to know 
"how maximal" 6*23 actually is. Of particular importance in this respect is the question 
whether the atmospheric neutrino mixing is exactly maximal because this would point to 
the presence of an underlying symmetry. This deviation will be known to order 10% after 
the next generation long-baseline experiments H3l . Comparative constraints are expected 
from future SK atmospheric neutrino data with statistics 50 times the current SK statistics 
|44j . Another way to state this is that we can identify if ^23 — 7r/4 is of order a/R or smaller 
with next generation long-baseline experiments. A deviation from 7r/4 even smaller (i.e., 
order R) can be achieved by dedicated next generation experiments HH HH] (see also 

m)- 

The effective mass in neutrinoless double beta is crucially dependent on the values of 
the mixing parameters 612 and 613. The current uncertainty on the values of sin^ ^12 is 
29% at the 3a level while sin^ 613 is completely unknown. The uncertainty on the value 
of sin^ 9i2 and the limit on sin^ 613 is however expected to reduce in the future. While the 
conventional beam experiments, MINOS, ICARUS and OPERA are expected to provide 
moderate improvement on sin^ 613 [42j, the Double-Chooz reactor experiment in France is 
expected to reduce the upper limit to sin^ 6^13 ^ 0.008 at the 3a level This upper limit 
could be improved further to sin^ ^13 ^ 0.0025 by the combination of the next generation 
beam experiments T2K j4S| and NOj/A |19|, as well as by the second generation reactor 
experiments ^Uj- It should be borne in mind that any one of these experiments could 
even measure a non-zero ^13, if the true value of ^13 happens to fall within their range of 
sensitivity. 

On the other hand, the improvement expected for sin^ 6*12 from the currently running 
experiments is small [HI |^ . The results from the third and the final Helium phase of 
the ongoing SNO experiment are expected to reduce the uncertainty on sin^ Ou to not 
more than ~ 21% and the KamLAND experiment is not expected to make any 

■^We define uncertainty of a quantity as tlie difference of the maximally and minimally allowed value 
divided by its sum and then multiplied by 100. 



6 



significant improvement on it ^HI^ESl- Doping the Super-Kamiokande (SK) detector 
with 0.1% of Gadohnium could improve our knowledge on the true value of sin^ 6*12 to 
18% uncertainty jHSj. The proposed/planned future experiments aiming to measure the 
very low energy pp flux coming from the sun are also expected to give a better measurement 
of the solar mixing angle jHSI- However, even with very small experimental uncertainty 
of only 1%, they are not expected to reduce the uncertainty on sin^6'i2 much better than 
~ 14% at 3(7 I56j. The only type of experiment that would provide extremely good 
measurement of the solar mixing angle is a long baseline reactor experiment with its baseline 
tuned to a Survival Probability minimum [SI], which is given by the condition Lmin — 
1.24 E/MeV eV^/Amj^. This type of experiment could be used to measure the mixing 
angle sin^ 6*12 down to ~ 6% at 3a jmEI]. 

The Dirac CP phase 6 has a faint chance of being measured in the next generation 
superbeam experiments jlHllin], provided the true value of 6*13 is not too small. However, 
performing an unambiguous measurement, and especially establishing a signal for CP 
violation, would probably need a beta beam facility or a neutrino factory ^58^. The two 
Majorana phases are measurable only in processes which violate lepton number. At present 
the only such process which seems to be viable experimentally is neutrinoless double beta 
decay jSH]. For detailed analyzes of how to extract information on Majorana phases from 
OuPP we refer to [inillHlliniEni- The bottom-line of this subject is that for given oscillation 
parameters the Majorana phases should not be too close to or 7r/2 and the uncertainty 
on both the experimental and theoretical side (i.e., the nuclear matrix elements) should be 
small fnil2ni- In addition, the prospects of determining the Majorana CP phases increase 
with increasing solar neutrino mixing angle. 

The quest for the limit on the absolute neutrino mass scale will witness attacks both via 
direct kinematical searches and cosmological measurements. The KATRIN |60| experiment, 
currently under construction in Germany, is scheduled to start taking data in 2008 and is 
sensitive to neutrino masses down to 0.2 eV. Further cosmological probes will test the sum 
of neutrino masses down to 0.1 eV within this decade Additional data sets and novel 
experiments can reduce this number by a factor of two [iil]. 

Finally, the last piece of information needed to construct completely the neutrino mass 
matrix is the ordering of the neutrino states - the sign of Amg^. If the sign was positive we 
would have a normal mass ordering, while for a negative sign we would have an inverted 
mass ordering. Typical approaches to identify sgn(Am|^) rely on using matter effects. 
Since the earth matter effect for modest baselines, and therefore modest matter densities, 
depends crucially on the size of the mixing angle 6*13, it becomes exceedingly difficult to 
study the mass hierarchy as 613 becomes small. If ^13 was large, close to its current limit, 
there could be a small chance of measuring the mass hierarchy using the synergies between 
T2K and NOi/A jH^j. However, the measurement would still not be very unambiguous 
and one would need either a beta beam facility or a neutrino factory for the hierarchy de- 
termination ^58j. Very large matter effects in the 1-3 channel are expected for supernova 
neutrinos. Therefore, a supernova neutrino signal could in principle be used to determine 
the sign of Arngj^ jHSl- Resonant matter effects in the 1-3 channel are also encountered by 
atmospheric neutrinos as they cross large baselines in their passage through the earth. This 
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can be exploited to probe the neutrino hierarchy both in water Cerenkov and large mag- 
netized iron calorimeter detectors jl3 IHl] • Recently some novel ways of probing the mass 
hierarchy requiring very precise measurements and using the "interference terms" between 
the different oscillation frequencies have been proposed |^ (see also [HE!)- However, it 
is understood that among all the neutrino parameters, the mass hierarchy determination 
is expected - together with the determination of the Majorana phases - to be the most 
challenging for the future experiments. In this respect the neutrinoless double beta decay 
is of some interest, since the question of distinguishing the neutrino mass hierarchy can be 
answered by Ou^p. This issue will be discussed in the following Section. 



3 Distinguishing the Neutrino Mass Schemes 

In this Section we present the phenomenology of neutrinoless double beta decay in terms 
of the oscillation parameters Arri^^, Amli, sin^ 6'i2 and sin^6'i3. In particular we will look 
how easy or difficult it would be to distinguish the different neutrino mass schemes if we 
have a signal or a significantly improved limit for Oi>(3(3. 

For the NH scheme, for mi -C m2 -C and therefore assuming that mi can be 
neglected, we have 



/ \NH 

(m) ^ 



sin^ Ou cos^ ^13 + ^Am^^ sin^ ^13 e^^^^"") . (13) 



The maximal value, {m)^^, is achieved when both terms add up, which corresponds to 
a — (3 = ±27r, ±7r or 0. The minimal value of (m) corresponds to a — (3 = ±vr/2, 37r/2, 57r/2 
or 77r/2, resulting in (partial) cancellation between the two terms in Eq. (I13|) . For non-zero 
values of mi between 10~^ and 10~^ eV there can be complete cancellation resulting in a 
vanishing effective mass. This will be subject of Sectional 

For the IH scheme, assuming that m3 ^ mi < m2 and neglecting m3, we have 



(m)'" ~ y I Am|i I cos^ ^13 VI - sin^ 2^12 sin^ a . (14) 

Here the maximal (minimal) value is obtained when a = 0, vr or 27r (a = 7r/2, 37r/2, 57r/2 
or 77r/2). 

Finally, for the QD mass spectrum 

{mf^ ~ mo I (cos' ^12 + sin' ^12 e''°) cos' ^13 + e^"^ sin' ^13] . (15) 

If the terms proportional to e'*" and e'*^ add up (i.e., when a and (3 take values of 0, vr 
or 27r) we have the maximal value of (m), which is then just rriQ. On the other hand, the 
minimal (m) is achieved when a and [3 take values of 7r/2, 37r/2, 571/2 or 77r/2. We will use 
this to set limits on mo in Sectional 

Hence, for a given mass scheme, the effective mass depends on Am|^, Am'i, sin' ^12, 
sin' ^13 and on the Majorana phases (to be precise, on one of them, or on a combination 



8 



thereof). In case of QD the absolute neutrino mass scale mo plays a decisive role as well. 
In the case of NH the smallest neutrino mass could be extremely important in deciding 
the degree of cancellation between the different terms, as discussed above. The lack of 
knowledge of most of those parameters means that we can not give definite predictions for 
the value of (m) for any of the three extreme neutrino mass schemes. We can however give 
a range of (m) for NH, IH and QD, namely: 

0.0 < (m) < 0.007 eV , NH (16) 

0.006 eV < (m) < 0.065 eV , IH (17) 
0.41 (0.07) eV < (m) < 2.3 (0.4) eV , QD (18) 

In order to obtain those values, we took the ranges of the oscillation parameters from 
Eqs. ()7p8p9piH) and inserted them in the expressions for the maximal and minimal values 
for (m) in the three extreme schemes under consideration. For the case of QD we used 
the limit from the direct kinematical search IHH] and in brackets mo < 0.4 eV. Note 
that with the 3cr values the NH and IH cases slightly overlap. However, these ranges are 
expected to reduce remarkably in the future from more precise measurements of the solar 
and atmospheric neutrino parameters and from either a signal for a non-zero value of 
(or from better upper limits on) sin^6'i3. Nevertheless, already at the present stage the 
possibility of distinguishing the mass spectra from each other opens up. In particular, the 
case of deciding between NH and IH is the most interesting one. Distinguishing the QD 
spectrum from NH or IH will be at most a consistency check since the common neutrino 
mass of or above 0.2 eV will be probed by either direct searches or cosmology. 

We present in Figure Q the lines of constant maximal and minimal (m) for the cases of 
NH and IH, respectively, in the sin^ 6'i2— sin^ ^^13 plane. We fixed the values of Am|^ = 0.002 
eV^ and Am^i = 8 x 10~^ eV^. The dotted (red) lines show the lines of constant minimum 
value of (m) for the IH scheme ((^)™n) and the solid (blue) lines give the lines of constant 
maximum value of (m) for the NH scheme ((^)max)- '^^^ Figure illustrates the well-known 
fact that (m)l^j^ has a strong sin^ 6*12 dependence, while (m)^|^ depends on both sin^ 6*12 and 
sin^6'i3. Also shown by the dot-dashed (black) line on the Figure are the current allowed 
values of sin^ 612 and sin^ 61^ at the 3ct level, obtained from a 2 parameter fit of the global 
oscillation data [10^^. As an example, one can see from the plot that if we know that the 
mass ordering is normal and assume that the smallest neutrino mass is negligible, values 
of the effective mass larger than roughly 0.0053 eV are incompatible with the currently 
allowed values of sin^ ^12 and sin^^i3. The same would be the case when we know that 
the mass ordering is inverted, assume that the smallest neutrino mass is negligible and the 
effective mass is larger (or smaller) than 0.023 (0.008) eV. 

We can see from Eqs. fjl3|) and ()14|) that the maximal and minimal values of (m) are 
crucially dependent on the values of the mass squared differences, especially on /S.m\^. 
Since we are talking about a future measurement of the Oz//5/3 signal, we take the expected 

■^The 3ct limit sin^ 6*13 < 0.044 given in the Introduction and Section El was obtained in a 1 parameter 
fit of the global neutrino data and is therefore different from the one in Figure ^ 
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futuristic uncertainties mentioned in Section El on Am|^ and Arn^i to plot the iso-range 
of and {fn)^^ in Figures El and O We can see clearly from the Figures that the 

uncertainty on the predicted value of (m) due to the uncertainty in the value of Am|]^ and 
Am^i would become very small in the next ten years. Therefore, in what follows, we will 
keep Amgi and Aml^ fixed at Aml^ = 0.002 eV^ and Am^i = 8 x 10~^ eV^. 

Any future positive signal for Oz//5/5 will be able to distinguish the IH scheme from 
the NH scheme if the measured value of (m) would be such that it could be explained by 
the predicted (m) for just one of the schemes and not the other. In other words, if the 
difference between the predicted values for (m) for the IH scheme and the NH scheme is 
larger than the error in the measured value of (m), then it would be possible (assuming 
that the smallest neutrino mass is negligible) to experimentally find the neutrino mass 
hierarchy from a measurement of Oz//?/3. We therefore show the difference in the predicted 
value of (m) for the two schemes, 

^{m) = - {m)ll , (19) 

in left-hand panel of Figure |31 as a function of sin^ 6*12 for 4 fixed values of sin^ ^13 and in 
right-hand panel of Figure (Has a function of sin^ ^13 for 4 fixed values of sin^ 612. We see 
that A(m) displays a strong dependence on sin^ 9i2, whereas the dependence on sin^ Oi^ is 
rather moderate. In fact, it holds that 

^ cos^ ^13 _ ^^^2 ^i + R + tan^ ^13) - tan^ ^13) , (20) 

1 + tan^ "12 

where R is the ratio of the solar and atmospheric Am^ defined in Eq. p2|) . Figure 01 and 
Eq. ()20p demonstrate the well-known fact that distinguishing the normal and inverted hi- 
erarchy is easier for smaller values of sin^ 612. The values of A(m), which are of the order 
0.01 eV, represent a value of the maximal experimental uncertainty which an experiment 
should have in order to be able to distinguish NH from IH. 

Up to here we neglected the complications which arise from the nuclear matrix elements 
(NME) involved in the calculation of the rate of Qvf3(3. There is neither a consensus in 
the literature on how large this uncertainty is, nor if there really is a large uncertainty 
at all [B7J. Let us therefore take a conservative point of view and analyze the situation 
as a function of the NME uncertainty. If there is indeed no problem with the NME, the 
statements given up to this point apply. 

Given the possibility of large uncertainties on the values of the nuclear matrix el- 
ements, it is plausible to ask if one could possibly extract the true hierarchy of neu- 
trino masses from a positive signal for Qvl3j3 in the future. To that end, we look for the 
difference between (m)™^ and {ra)^_^ after including the uncertainties coming from the 
nuclear matrix elements. A way to perform such an analysis has been developed in |16j . 
and here we follow a very similar approach: in order to parametrize the uncertainty com- 
ing from the nuclear matrix elements, we define the {^v(3(3 decay rate F measured in any 
experiment such that 

Vf = x(m), (21) 
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where x contains all other factors involved in the Oz//?/3 decay, including the NME. Thus 
for any measured T the measured value of (m) would be 

("^)measurcd = VT / X ■ (22) 

If y is the smallest possible value for x, and if we assume that the uncertainty on x originates 
solely from the uncertainty on the NME and is a factor of z (with z > 1), then the range 
of X is given hj y < x < zy. Therefore, the uncertainty on the measured value of (m) that 
we have due to the uncertainty on the NME is given by 

z — 1 

^("^)measured = • (23) 

In order to be able to experimentally distinguish IH from NH with the help of a Ou^f] 
measurement, we must have {fn)^^^^ > (m)^^ after including the uncertainty from the 
NME. Thus we need the condition 

- ^ Hll > (24) 

to ascertain the neutrino mass hierarchy. In the left and right panel of Figure El we plot the 
difference — z {m)^^ predicted as a function of sin^ 6^12 and sin^6'i3, respectively, 

including the uncertainty coming from the NME. We show the plots assuming a range of 
values for the NME uncertainty z. As a function of sin^ 6^12 the quantity {fn)™^ — z {m)^^ 
shows a similar behavior as for no uncertainty in the NME. The dependence on sin^ ^13 is 
rather weak for small z, but becomes larger with increasing z. From Figure El one can read 
of the maximal experimental uncertainty for a given z, sin^ ^13 and sin^ 9i2 for which one 
could still be able to distinguish the mass hierarchy. For instance, taking sin^ ^13 close to 
its current hmit and assuming z = 2 and sin^ 6^12 = 0.3, we have — z (fn)^^ ~ 0.01 

eV and the experimental uncertainty should not be larger than this value. For zero 6*13 
this value becomes larger by roughly 40%. 



Up to now we considered the case of the smallest mass being zero, which would 
correspond to the neutrino mass matrix having a zero determinant, clearly a rather limited 
and special case. For a given set of oscillation parameters and the nuclear matrix elements 
uncertainty, it would be possible to distinguish between the NH and IH in a Ouf^f] experi- 
ment only if is less than a certain value. It is interesting to examine up to which values 
of the lightest neutrino mass one can distinguish NH from IH. Toward this exercise, we 
plot in Figure (HI the difference between the minimal value of (m) for IH and the maximal 
value of (m) for NH, {m)™^ — {m)^^, as a function of the smallest neutrino mass msm- 
The NME uncertainty is also included in the analysis. For each of the 4 panels we assumed 
a certain "true" value of sin^ Ou with an optimistic uncertainty of 6% jUlEZI- For sin^ 613 
we assumed that no positive signal is observed in any of the forthcoming experiments in 
the next ten years. Thus we allowed it to vary within and 0.0025 42 . We see from the 
Figure that for sin^ 6*12 = 0.3 for example, when z = 1.5 {z = 1, viz. no NME uncertainty) 
and A(m) = 0.01 eV, we can distinguish NH from IH if msm is smaller than 0.002 (0.004) 
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eV. Of course, the value of A(m) obtained like this is again an upper limit on the ex- 
perimental uncertainty which would still allow distinguishing NH from IH for a non-zero 
smallest mass. If sin^6'i2 = 0.26 the situation looks more promising. For the same value 
of A(m) = 0.01 eV we can distinguish NH from IH if rngm is smaller than 0.005 (0.01) eV 
ii z = 1.5 {z = 1). Other cases are easily read off Figure IHl 

For a larger uncertainty in the value of sin^ 9i2 and/or for a larger true value of sin^ ^13, 
the chances of distinguishing the IH from NH become worse and in addition only work for 
very small values of m^m. 

We see however that in principle values of rngm up to ~ 0.01 eV are allowed in order to 
distinguish the normal from the inverted neutrino mass hierarchy. 



4 Limit on the Neutrino Mass 

Another interesting aspect of neutrinoless double beta decay is the possibility to set a limit 
on the absolute scale of the neutrino mass. The current upper value for the effective mass 
is somewhere between 0.3 and 1 eV, where this range of course has its origin in the NME 
uncertainty. The indicated values correspond to the QD mass spectrum, on which we wish 
to focus in this Section. As indicated in Sectional for a given common mass scale mo, the 
lowest possible effective mass can be written as 



/ \QD /irr |2 irr |2 \tt \2\ " taU^ ^^12 — 2 | f/e3 P 



Hence, having a limit on the effective mass at hand, we can translate it into a limit on the 
neutrino mass. We write the experimental limit as 

(m)-P = z {mTZ , (26) 

where (m)|^? is defined as the limit on (m) obtained by using the largest available NME 
and z > 1 encodes again the NME uncertainty. Then, the limit on the neutrino mass reads 

We have introduced here a function f {612, (^13) in this expression, which separates the 
information available from neutrino oscillation experiments from the information coming 
from Oz//5/5 and mo. We think that this formulation might be helpful in understanding how 
neutrinoless double beta decay and the absolute neutrino mass scale are related. Currently 
the uncertainty on f{9i2,0i3) is around 50%, 1.9 < f{9i2,9i3) < 5.6. It is expected 
to reduce to ~ 21%(~ 9) % at 3cr if a low energy pp solar neutrino experiment (reactor 
experiment at the survival probability minimum) would be built. The uncertainty depends 
only little on the value of ^13. From the current limit on the effective mass, (m) < 0.35 z 
eV, with the accepted value of 2 ~ 3 for ™Ge (see ^J), and our current knowledge of 
f {012, 613), we can set a limit on mo of 5.6 eV, clearly weaker than the limit from tritium 
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experiments. Only for close to vanishing NME uncertainty (i.e., z = 1) we can reach values 
of mo ^ 2 eV which are then comparable to ones from direct searches. 

In Figure 13 we show the iso-contours of mo predicted for the QD mass spectrum (cf. 
Eq. (fTTjl ) in the sin^ 9i2 — sin^ 6'i3 plane. For each of the iso-contours we have assumed an 
illustrative value of z (m)^;"^ = 0.1 eV. Also shown is the current 3a allowed area of sin^ 6^12 
and sin^ 613. From the Figure we can see that for a measured value of z (m)^j'^ = 0.1 eV, 
the constraint on the common mass scale for a QD spectrum would be: 0.2 eV ^ mo ^ 0.6 
eV. For any other measured value of z (m)^? corresponding to a QD mass scheme, we 
can obtain the corresponding constraints on mo by scaling the above limit suitably, i.e., 
for z (m)l^? = 0.2 eV we would have 0.4 eV ^ mo ^ 1.2 eV. The limits on mo given 
above are of course with our existing knowledge about 6*12 and ^^13. With reduction of the 
range of allowed values of 6^12 and ^13, especially 6*12, the constraints on mo are expected 
to reduce substantially, as discussed above. For instance, if /(^i2,^i3) was known with an 
uncertainty of 20%, say 2.7 < /(6'i2, 6*13) < 4.0, then we would have for z (m)^? =0.1 (0.2) 
eV that 0.3 eV ^ mo ^ 0.4 eV (0.6 eV ^ mo ^ 0.8 eV). Of course, if we have no signal for 
OuPP, but just an upper limit on z {m)^-^^, we have no longer an allowed range on mo, but 
an upper limit corresponding to the largest value in the range. From the examples given 
above, one can note that for the QD mass spectrum, a measurement or a better constraint 
on (m) will lead to a stronger limit on the absolute neutrino mass scale compared to the 
current limit from direct kinematical searches. 

5 Implications of a Vanishing effective Mass 

It is well-known that not only a measurement but also a non-measurement of neutrinoless 
double beta decay has significant influence on our knowledge of the unknown neutrino 
parameters. In this Section we wish to discuss a rarely studied subject, namely the impact 
that a very small upper limit on (m) would have^. An exactly vanishing (m) would 
correspond to a texture zero in the neutrino mass matrix in the charged lepton basis - 
certainly an interesting feature. However, to prove/observe an exactly vanishing effective 
mass is a formidable task. 

In case of an inverted hierarchy we know that there is a lower limit of (m) which is 
— when using current 3a values of the oscillation parameters — given by roughly 0.006 
eV, see Eq. (fTTjl . Let us assume that neutrinos are Majorana particles, have an inverted 
ordering and that the effective mass takes has an upper limit of order 0.01 eV. This is 
a situation which might arise if we know from some other independent experiment that 
sgn(Am^) = —1 and the next generation 0i//9/9-experiments do not find a signal corre- 
sponding to a non-vanishing (m) , but give an upper limit on (m) still above the theoretical 
limit (m)™j^. Then we can infer from Eq. (fT^ the values of sin^ 9i2 and the Majorana phase 
a which are still compatible with the data. In Fig. |H1 we display the result for an expe- 
rimental upper limit on (m) of 0.04, 0.03, 0.02, 0.01 eV, taking m3 = 0, ^13 = and 
Amg^ = 0.002 eV^. We checked that the results are rather stable when we depart from 

^For a related earlier analysis, see |68j . 
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these values within their current uncertainty. The current 3a values of sin^ 6*12 are also 
indicated in the Figure. Allowed is the area to the right of the respective curves. For 
instance, if (m) < 0.02 eV and sin^ 612 = 0.3 then a has to lie between 1.2 and 1.9, or 
a ~ 7r/2 ± 0.4. Alternatively, if (m) < 0.02 eV and sin^ 6^12 < 0.28, then the IH case is 
ruled out. 



Figure ini shows a similar analysis for the case of NH with a smallest mass mi = 0. In our 
parametrization, as can be seen from Eq. (fTSjl . it is the combination j3 — a which governs the 
destructive interference which leads to very small or zero (m). We took very small upper 
limits on (m) of 0.004, 0.003, 0.002 and 0.001 eV (hence this discussion will not become 
realistic within the next 10 years) and chose mi = 0, sin^ Oi^ = 0.04 Amgi = 8 x 10~^ 
eV^ and Am|i = 0.002 eV^. Note that for very small values of 613, corresponding to 
sin^6'i3 <^ R sin^ 6*12 ~ 0.01, the dependence on this combination of phases drops. In this 
part of the Figure it is the left of the respective curve which is allowed. For instance, for 
(m) ^ 0.002 eV and a rather large value of sin^ 9^ = 0.4 the phases have to be very close 
to 7r/2, namely /3 - a ~ 7r/2 ± 0.2. Alternatively, if (m) < 0.001 eV and sin^ 6112 > 0.33, 
then the NH case is ruled out. 

As mentioned earlier in Section El we note here that in the case of NH, (m) could have 
a sizable dependence on the smallest neutrino mass state mi. We plot therefore in Fig. ITOl 
the same as in Fig. IHlbut for a smallest neutrino mass of mi = 0.005 eV. We took again 
sin^6'i3 = 0.04, Am^i = 8 x 10"^ eV^ and Am|i = 0.002 eV^ and chose the same four 
upper limits on (m) as above. With a non-vanishing mi Eq. ()13p is modified to 



(28) 



(m) ^ mi cos^ 6112 cos^ 6113 + yjml + Am| sin^ 6^12 cos^ 6^13 e^*" + \j iAm\ sin^ 6113 e^*'^ 

and thus it is no longer (3 — a but both a and (3 which play a role. In Figure ^1 we 
fixed a = tt/2 which gives a negative sign for the second term in Eq. (j^H|l . We indicate 
the allowed (disallowed) values of the parameter regions in the Figure with A (D). If 
(m) < 0.004 eV and (m) < 0.003 eV the disallowed areas are already outside the currently 
allowed 3(T region of sin^6'i2 for a fixed at 7r/2. To show the dependence of our results on 
a, we plot in Figure ^2 the allowed areas in the a-j3 plane for four different fixed values of 
sin^ ^12, keeping sin^ 6^13 fixed at 0.04 and for mi = 0.005 eV. Allowed are the areas within 
the respective curves. 

An interesting aspect is when ^13 vanishes but mi is non-zero. Then the effective mass 
is a function of the phase a alone: 



(m) 



mi cos^ 6^12 + ^/m? + Am% sin^ 6^12 e^*° 



(29) 



In Figure ^1 we show the areas in the sin^ 6'12-a plane that would be allowed if we had an 
upper limit on (m) of 0.001 eV (black line), 0.002 eV (red line), 0.003 eV (green line) and 
0.004 eV (blue line). The areas outside the closed curves would be disallowed. We can see 
that for a limit of (m) =0.001 eV, all values of a could be possible while for (m) =0.004 eV, 
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only the range 7r/2 ± 0.67 is allowed. This allowed range of a is not expected to improve 
with the reduction of the uncertainty on sin^ 6^12 since its 6*12 dependence, as can be seen 
from the Figure, is very weak. 

For values of the smallest mass of m\ ~ 0.005 eV and sin^ dyi ~ 0.3 both terms in Eq. 
are of roughly the same magnitude (i.e., m\ ~ m2) and we could write 

(m) ~ mNH \/l — sin^ 2Qyi sin^ a , where = + Am| ~ 0.01 eV. 

Note the similarity of this equation with Eq. ()14|) in case of IH. 

To sum up, in certain cases it might be possible to significantly constrain the allowed 
values of the Majorana phases. Moreover, the allowed values of the solar neutrino oscillation 
parameter Q\2 can be comparable to its current known range. 

6 Neutrinoless Double Beta Decay, Future Neutrino 
Oscillation Data and the Identification of the Neu- 
trino Mass Matrix 

In this Section we wish to give a summary of typical neutrino mass matrices available in 
the vast literature. Predictions for and correlations between the neutrino observables are 
implied by each of the candidates and can be used to distinguish them. In particular, 
knowledge of the neutrino mass spectrum, i?, ?7e3, 6*23 and of course the effective mass are 
then very helpful for identifying the correct mass matrix. A complete study of all possibil- 
ities is (and maybe can) not performed, however, we feel that the majority of the models 
will produce at the end of the day a mass matrix which will more or less correspond to 
one of the examples given here. For instance, there are many symmetries which eventually 
lead to a /ir symmetric mass matrix [69j, see for instance Refs. [7D] for some examples. For 
more models and details we refer to the excellent reviews available [71]. We should remark 
that for the given examples, unless otherwise stated, the charged lepton mass matrix nii is 
diagonal, i.e., for the matrix \Ji diagonalizing m^mj holds that Ui = 1, and therefore the 
PMNS matrix is just U = Uj = U^, where Ui, diagonalizes m^. Further note that most 
correlations focus on the quantities (m), R, 623, If/esI or 6^12. Interplay with the parameters 
describing CP violation is rarely studied. 

Let us focus first on the normal hierarchy, for which there is a very successful and 
well-known texture, namely [72] 



/ ae 



:2 be de \ 



my = mo 



(30) 



V 



9 I 
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Here a, b, d, e, /, g are complex parameters of order one and e is a real and small parameter 
typically of order of the Cabibbo angle A ~ 0.22. Such a mass matrix preserves at lead- 
ing order the lepton charge Lg and can be obtained already by simple models based on 
U{1) charges or through sequential right-handed neutrino dominance [Z3I- A correlation 
resulting from Eq. (jHUjl is 



(m) = ci ^JAml \Ue3\^ with |f/,3| = Ca , (31) 

where Ci^2 are functions of the order one parameters and naturally also of order one. With 
the parameters in the fir block unspecified, the atmospheric neutrino mixing angle typically 
deviates sizable from its maximal value 7r/4, i.e., 6'23 = 7r/4 — C3 \/R. 

The discussion is also possible in the context of symmetric mass matrices [^11 EO] . 
Since their predictions include maximal 623 and zero ^13, one is interested in breaking 
schemes of the symmetry. Breaking the symmetry in the e sector of leads to the same 
correlation [73] as in Eq. (|3ip but atmospheric neutrino mixing is much closer to 7r/4. On 
the other hand, breaking the /xr symmetry in the fir sector of rrii, leads to [7i] 

{m) ci ^fX^lUesl , (32) 

where {Ue^l is now very small, |t/e3| = C2-R, and the deviation from maximal atmospheric 
neutrino mixing is again sizable, i.e., of order \/R. Another interesting breaking scenario 
of fiT symmetric mass matrices is when the e/z element is not identical to, but the complex 
conjugate of the er element [7^]. Then one finds that Ue3 7^ but Re f/gs = 0, i.e., there is 
maximal CP violation with S = ±7r/2. A special case of fir symmetric mass matrices, for 
which similar statements will apply, is given when tan^ 612 = 0.5, which then corresponds 
to the so-called tri-bimaximal mixing scheme [70] . "Bimaximal" scenarios [77] with their 
prediction tan^ 9^ = 1 typically require contributions from the charged leptons, see below. 

Perturbing a zeroth order mass matrix that leads to zero f/es and maximal 6'23 with 
a random matrix containing small entries of the same order e will lead to a semi- 
hierarchical mass spectrum (i.e., m2 > mi and not m2 ^ mi) with a somewhat larger (m) 
and also to sizable deviations both from Ues = and 6*23 = 7r/4 f7Hj . 

Other candidates studied frequently in the literature are minimal models in the sense 
of having zeros in mi, [70] or a minimal number of parameters in a see-saw context [80j . 

Table Q summarizes our collection of the typical mass matrices for the normal mass 
hierarchy. Once we know that the neutrino mass spectrum is hierarchical, an information 
that could be provided by a negative search for OujSP in future experiments and knowing 
that sgn(Am|^) = +1, we can sort out the correct mass matrix when we have precise 
information on the oscillation parameters. 

In case of the inverted hierarchy, summarized in Table El the most stable candidate 
"theory" corresponds to a mass matrix conserving the flavor charge — — Lr |81j . 
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Such a matrix is given by 



/ cos 6* sin6' \ 



m„ = mo 







V 





/ 



(33) 



predicting one massless neutrino, zero di^, zero (m), maximal so/ar neutrino mixing with 
Attt-q = and 623 = 9. Barring extreme fine-tuning, it is impossible to perturb the 
structure of Eq. (jHHj) such that the measured values of 6*12 and Am^ are both accommodated 
in accordance with the data. Typically, after breaking it will hold that tan^ 612 — 1 ± i?/2 
(sin^ 6'i2 — I ± -R/4), whereas experimentally sin^ 612 — \ — VR is required. An exception 
is a see-saw model in which the perturbations at high energy have the same order of 
magnitude as the terms allowed by the symmetry, see jH2- In that example, Ue^ and the 
smallest mass state remain zero. Another way out is to take corrections from the charged 
lepton sector into account |83,. The matrix Ui, which diagonalizes mim\, is multiplied 
from the left to f/jy, which is the matrix diagonalizing my from Eq. (jHHj) . The latter has 
to be perturbed in order to generate a non-zero Ami,. Given the observation that the 
deviation from maximal solar neutrino mixing is determined by the Cabibbo angle, i.e., 
7r/4 - ~ A ^ 0.22 ~ Oc il, one assumes a CKM-like structure of jHS]^- The 
typical result of the described Ansatz lies in a correlation between the solar neutrino 
mixing and ?7e3, namely sin^ 6*12 — \ — cos0 cot 6*23 \ Uez\ , where is a phase appearing 
in rriu, which can be given by the Dirac phases as measurable in oscillation experiments. 
This depends however on the breaking in niy. Nevertheless, both If/esI and tan^ ^12 are 
expected to be close to their current upper limits. Assuming minimal breaking in m,^ 
[SB] one can show that there are remarkable correlations between the observables, such 
as (m) ~ \/ Am^ | cos 2^12 + Ai Jcp/ sin^ 6'23|, where Jcp is the Jarlskog invariant for CP 
violating effects in neutrino oscillations, which is proportional to |f/e3|/4. This represents 
one of the few examples where the CP phases are part of the predicted correlations, 
depending however on the details of the model. The correlation tan^ ^12 ~ 1 — 4 |f/e3| is 
however independent of the breaking in and relies only on bi-large jSH]- Hence, the 
two discussed possibilities incorporating L^ — L^ — predict values for (m) of similar size 
but predict either zero or large Ue^- 

There is another zeroth order scheme for the inverted hierarchy, namely 



/ 1 \ 



m,. 



mo 



\ 



1/2 1/2 
• 1/2/ 



(34) 



Zero f/e3 and maximal ^23 are predicted and the two leading mass states have the same 
sign. One can perturb this structure [TS] and the result is that both Ues and ^23 — 7r/4 are 

^Note that the so-called Quark-Lepton-Complementarity, i.e., the exact relation 6*12 + 6c = Tr/4: is a 
special case of this procedure. 
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very small quantities, but the effective mass is larger than in case of matrices based on the 
conservation of — — L^-. 

The dependence of (m) on Ues goes from none to sizable to little in the three cases 
discussed, cf. with Table El So, if we know that the mass ordering is inverted, i.e., 
sgn(Am3]^) = — 1 and the effective mass is of order a/ Am^, then we basically only have to 
distinguish between Lg — — from Eq. (jH^ and the matrix given in Eq. (jS^)- This 
can be done by knowing how to close (m) is to \/ Am\. More precision data on 1/^3 and 
$23 will fix the details of the model. 

Now we turn to quasi-degenerate neutrinos. Usually, in such scenarios both the normal 
and the inverted mass ordering can be accommodated. Typical examples are given in Table 
El There are matrices with two zero entries which are compatible with quasi-degenerate 
neutrinos and which predict an interesting interplay of observables [70]. Anarchical mass 
matrices ^] do in general not generate extreme mixing angles and therefore typically 
result in large deviations from zero 11^3 and maximal 623, i.e., of order R. 

Possibilities to obtain QD neutrinos are models based on flavor democracy jHZI, of which 
we quote in Table (21 one interesting recent and typical example |HH] • In such models both 
f/j, and Ui are required in order to reproduce the neutrino data, resulting in a dependence 
of the neutrino observables on the charged lepton masses. This makes the predictions 
also a function of the democracy breaking scenario, therefore somewhat model-dependent. 
However, the smallness of \Ue3\ in such models can be attributed to the small ratios of the 
charged lepton masses. 

One can also upgrade a hierarchical mass spectrum (e.g. resulting from sequential dom- 
inance) to a quasi-degenerate one by adding a term proportional to the unit matrix (e.g. 
connected to a SO (3) symmetry) to it [11]. The neutrino mixing angles display approxi- 
mately the same behavior as in the NH case (i.e., as from Eq. ((201)), however, the effective 
mass is close to the common mass scale, i.e., there is little cancellation in (m). Moreover, 
the larger mo the smaller the phases [5U] . 

Some attention has recently been caught by a model based on the discrete symmetry A4 
|90j . Applying the most general SUSY threshold corrections to the resulting mass matrix 
leads to predictions such as very large atmospheric mixing, purely imaginary Ue3 with 
an absolute value of order \/R and (m) = mo. Similarly, one can use a simple Abelian 
symmetry corresponding to the conservation of the flavor symmetry — Lr [HII and 
perturb the resulting mass matrix'' 



= mo 



/ cos9 \ 

sin^ 
V • • / 



(35) 



to which the A4 model corresponds when 6 = tt/A. The deviations from zero f/e3 and 
maximal 623 are small and inverse proportional to mp. Note that interestingly most of the 

^Note that this matrix is a special case of /ir symmetry. 



18 



models presented here have (m) ~ mo, i.e., the mass parameters hopefully measured in 
direct laboratory searches and in neutrinoless double beta decay experiments should be 
almost identical. 

Another zeroth order scheme for quasi-degenerate neutrinos corresponds to a matrix 
of the form |92J 



/ -I/V2 1/v^ \ 



rrii, = mo 



V 



+ (l+r/)/2 

(l+^)/2/ 



(36) 



Exact bimaximal mixing is predicted and breaking of the matrix is also required to generate 
splittings between the mass states. Similar statements as for the matrix Eq. (j33|l hold, i.e., 
without extreme fine-tuning it is impossible to generate deviations from maximal solar 
neutrino mixing sizable enough not to be in conflict with the data. Hence, the contribution 
from the charged lepton sector are required, which will lead to similar correlations as for 
Lg — — Lt- discussed above, such as tan^ 6'i2 — 1 — 4 \Ue3\- The correlation between the 
other observables depends strongly on the breaking. Since however the ee element is zero 
in Eq. (jSE)), the effective mass is expected to be small, i.e., (m) ~ mo cos 26*12. 

Hence, the identification of the mass matrix in case of a QD spectrum can be performed 
when information also from direct kinematical or from cosmological measurements is avail- 
able. With the expected future precision data the identification of the mass matrix should 
also be possible. 



To sum up this Section, a complete determination of the all neutrino observables and 
consequent identification of correlations will be able to discriminate between the various 
possible models. At least some of the possible candidates will be ruled out. Note finally 
that the values of the two Majorana phases, whose determination turned out to be an 
extremely challenging task, is not really required in order to identify the mass matrix. 



7 Conclusions 

We have analyzed some aspects of the connection between neutrinoless double beta decay, 
the neutrino mass scale and neutrino oscillation parameters. In particular, we concentrated 
on the question whether the expected future precision on the oscillation parameters sim- 
plifies certain important consequences of a measurement of, or improved upper limit on, 
Ou^p. We first summarized the current situation of the determination of the nine physical 
parameters of the neutrino mass matrix and the prospects of improving our knowledge 
about them with future experiments. Then we analyzed how Oi>(3f3 can help in distin- 
guishing the normal mass ordering from the inverted one. We included the nuclear matrix 
element uncertainty in the analysis and pointed attention to the fact that distinguishing 
NH from IH depends on the value of the smallest neutrino mass rrism- We analyzed this 
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point and found that in principle values of rrism = 0.01 eV are allowed, where of course the 
issue of the nuclear matrix elements can significantly spoil this possibility. 

Then we investigated inasmuch Ou^jS can be used to set a limit on the neutrino mass 
scale mo in case of a quasi-degenerate spectrum. We argued how the information from 
oscillation data and from Ou/SP can be separated. Current limits on mo are weaker than 
direct ones but can be significantly improved with future Oi//3/3-experiments and more 
precise knowledge of the oscillation parameters sin^6'i2 and sin^^i3. 

Next we studied the implications of a very small, or even vanishing effective mass. 
Still insisting that neutrinos are Majorana particles, one can put in this case interesting 
constraints on parameters, in particular on sin^ and the Majorana phases. 

Finally, we tried to perform a scan through the literature and identified neutrino mass 
matrices, which tend to arise frequently in many different models. We listed their predic- 
tions and correlations and argued how future precision date of the oscillation parameters 
and information from neutrinoless double beta decay can help to identify the true neutrino 
mass matrix or at least to sort out many unsuccessful ones. 
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Table 1: Typical and popular neutrino mass matrices predicting the normal mass hier- 
archy. Given are some comments about their origin and order of magnitude predictions 
for correlations between the observables. The parameter e is usually small, i.e., mostly of 
order A ^ 0.22 and the parameters a, 6, d, e, /, ci,2,3 are of order one. 
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Table 2: Same as Table ^ for the inverted mass hierarchy. We have Jcp = 
1/8 sin 2012 sin 2023 sin20i3 cos0i3 sin 5 with the neutrino oscillation CP phase 6. The 
correlations of the last row depend on the breaking of the mass matrix. For more compli- 
cated breaking, the phase 6 will be some combination of phases. 
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Table 3: Same as Tabled for the quasi-degenerate mass spectrum. The quantity measur- 
able in direct mass measurements is denoted mo. 




Figure 1: Lines of constant maximum value of (m) (in eV) for the NH scheme ((m) 
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NH 



0.1 



Am^3j=2.0xl0 eV" 
Am^3j = l. 91x10"^ eV^ 
Am^,, =2.09x10"^ eV^ 



^a=" 0.05 



O 



♦ • ^ 

♦ . * 

♦ • 4- 
% : t 

♦ ■• I 
% : » 

♦ - 

♦ • * 

♦ • ^ 

$ if- 



♦ • ^ 

♦ - 

♦ • ^ 

♦ ■ 



|H0.025 



t • |0.03 
S » •- * 

|0.035|":| 



♦ - « 

♦ - i*: 



- * 
• i 



♦ -^^^ 

♦ •* 
♦ 



|-* 0.02 

■ + 

lit 



/ 



U O.Ol 



|:t 0.015 



0.1 



0.2 



0.3 

sm 0,, 



0.4 



5 0.005 

a 



0.5 



Figure 2: Lines of constant minimal value of (m) (in eV) in the IH case for 
different Amoi . 
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Figure 3: Lines of constant maximal value of (m) (in eV) in the NH case (('^)max) 
different Amh . 




sin sin 6^3 

Figure 4: The difference between the minimal value of (m) for IH and the maximal value 
of (m) for NH, (^«)™^ — {raf^^, as a function of sin^6'i2 (left-hand panel) and sin^ 6^13 
(right-hand panel). 
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Figure 5: The difference between the minimal value of (m) for IH and the maximal value 
of (m) for NH times the NME uncertainty, — z (^)n^ax' ^ function of sin^^is for 

different z and sin^ = 0.28 (left-hand panel) and sin^^i3 (right-hand panel). 
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Figure 6: The difference between the minimal value of (m) for IH and the maximal value 
of (m) for NH, (m)™^^ — ^^jnf^y^-, as a function of the smallest neutrino mass m^^. The 
4 panels are for 4 assumed "true" values of sin^6'i2. The values of sin^6'i2 are allowed to 
vary within 6% of the assumed true value, while sin^ ^13 is allowed to vary between and 
0.0025. 




Figure 7: 
spectrum. 



The lines of constant uiq in the sin^ 9i2 — sin^ ^13 plane, predicted for a QD mass 
Also shown is the current 3a allowed region in sin^ 9^ and sin^ 613. 
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Figure 8: Allowed values of sin^ 9i2 and the Majorana phase a for a small upper limit on 
(m) (in eV) in case of IH. The dotted lines show the current 3a limit on sin^6'i2. The 
allowed areas are to the right of the respective curves. 
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Figure 9: Allowed values of sin^ and the Majorana phase combination P — a for a small 
upper limit on (m) (in eV) in case of NH. The dotted lines show the current 3a limit on 
sin^6'i2. The allowed areas are to the left of the respective curves. 
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Figure 10: Allowed values of sin^ 9i2 and the Majorana phase /3 for a smallest neutrino 
mass of mi = 0.005 eV in case of NH for 4 different assumed upper limits on (m) and 
a = 7t/2. The dotted lines show the current 3a limit on sin^ 6*12. The allowed (disallowed) 
areas are marked with A (D). 
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Figure 11: Allowed values of the Majorana phases a and (3 for a smallest neutrino mass of 
mi = 0.005 eV in case of NH for 4 different assumed upper limits on (m). The 4 panels 
correspond to 4 different true values of sin^ 612 shown on the Figure. Allowed are the areas 
within the respective curves. 
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Figure 12: Allowed values of the Majorana phase a and sin^ 6*12 for NH for sin^ 6*13 = 
and with a non-zero smallest neutrino mass of mi = 0.005 eV in case of NH for 4 different 
assumed upper limits on (m). The dotted lines show the current 3cr limit on sin^6'i2. 
Allowed are the areas within the respective curves. 



